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ABSTRACT 



The radiation impedance of a finite array of finite cylinders is 
determined by a boundary value problem technique. The final equations 
for the radiation impedance are stated without restriction on oylinder 
length or radius, oylindor spacing, or on the amplitude and phase of 
the radial velooity of the cylinders. The derivation is made for the 
case of each cylinder in tho array having a uniform radial velocity over 
a soctor of the oylinder, but tho method of solution is not limited to 
such a velocity distribution. The results of a numerical calculation of 
tho effect of electrical beam steering on the radiation impedance of a 
specific array are shown to illustrate such effect and the uso of the 
derived equation for radiation impedance. 

The writer wishes to express his appreciation for the assistance and 
onoouragoment given him by I.Ir. John l'liokman of the U. S. Ilavy Electronics 
Laboratory and Professor II. liedwin of the U, S. Naval Postgraduate School 
in this investigation. 
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Wave velocity 
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Instantaneous acoustic pressure 
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Volume density 
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The notation associated with the Bessel Functions used in this paper 
is taken from N. W. McLaohlan, Bessel Functions for Engineers, Oxford 
at the Clarendon Press, 1955* 



Jn<*> 

V*) 

w (D,(2) 



i 
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Bessel function of the first kind 

Bessel function of the second kind, as defined by Weber 
Modified Bessel function of the second kind 
Bessel function of the third kind# 

H (l h) = J (z) -f jY (z) 
n n n 

( 2 ) 

( z ) "* ^n( z ) O-^n^ 2 ) 

The positive square root of a real, positive quantity. 
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1# Introduction 



The radiation impedance of a finite array of finite cylinders in an 
infinitely long, rigid, cylindrical baffle has been investigated by D. H# 
Robey# ^ Robey considered the problem of uniformly pulsating cylinders 
and his final conclusions are based on a cylindrical array of radius 
negligible in comparison to the wavelength of the radiated sound# It is 
the intent of this paper to investigate the radiation impedance of such a 
finite array of finite cylinders with fewer restrictions imposed on the 
radial velocity distribution and size of the array# It is also the intent 
of this paper to apply the results of such an investigation to the study 
of the effects of electrical beam steering on the radiation impedance of 
a finite array of finite cylinders# 

The radial velocity distribution on each pulsating cylinder will be 
specified as a boundary condition to the wave equation which is to be 
solved# However this is not enough to completely specify a unique solu- 
tion to the wave equation# Therefore it is proposed that as the radial 
distance from the cylindrical radiator increases, vanishing boundary 
values for pressure and velocity aro required# This condition again 
proves too weak to determine a unique solution# Therefore a more stringent 
boundary condition in addition to the condition of vanishing boundary 
values is required, ie. n The Sommerfeld Condition of Radiation#* 1 That is, 
"The energy which is radiated from the souroes must scatter to infinity; 

^D# H. Robey, On the Radiation Impedanoe of an Array of Finite Cylinders, 
J# Acoust# Soo# Am#, 27, pp706-711, 1955 

^A# Sommerfeld, Partial Differential Equations in Physics, Academic 
Press Inc., 1949, p 189 
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no energy may be radiated from infinity into tho prescribed singularities 
of the field (piano waves are excluded since for them even the condition 
u - 0 fails to hold at infinity). 11 

The wave equation which is to be satisfied is stated in cylindrical 
coordinates : 

5r 2 r Sr r 2 5 e 2 d z 2 c 2 dt 2 
The scalar velocity potential function is understood to be ^(r,e,z,t). 

The z axis of the cylindrical system corresponds to tho longitudinal axis 
of the cylindrical array (see figure l). 

The driving forces which are to be imposed on the elements of the 
cylindrical array will vary harmonically with time. Therefore for the 
results of this paper to conform to standard electrical engineering 
practice of assigning positive reactance to inductance, a time dependence 
of the form will be assumed. That is j^(r, e,z,t) - 0(r,e,z) 

After substitution for ^(r,ez,t), tho wave equation is reduced to the 
following: 

4 JL - 4 £& 4 —1 4 £$ z 0 

> 2 r or 2 -.2 *2 2 

dr r de dz c 

A formal solution to the above equation is discussed in Appendix 1. 

The radiation impedance of an array of finite cylinders will now be 
determined. 
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2 0 Derivation 



The first array to be investigated is the least restricted in 
boundary conditions of the arrays to be consideredo The array is to be 
constructed as follows (soe figure l): 

1. There are m active cylinders of length = 2^*, radius = a , 
and center to center spacing = 2b* 

2 • The cylinders are mounted coaxially in an infinitely long, 
rigid, cylindrical baffle of radius = a# 

The pulsation of any active cylinder is to be described as follows: 

1* On the given cylinder, the radial velocity is constant 
over the length of the cylinder for a given aximuth angle* 

2. On all cylinders, the radial velocity is uniform between azimuth 
coordinates - <X to */ oC and zero elsewhere* 

Therefore the radial velocity of the i^* 1 cylinder may be described by the 

following Fourier Series: 

oo 



h - 




(2 o 0l) 



Where A s s n sin nc < 
n 


> e n 


= 1» 


o 

D 

£ 


ncc 


6 n 


= 2, 


n “ 1,2,3, •••• 



and - the phase of the pulsations. 

The scalar velocity potential 0(r, e, z,t ) must satisfy the following 
wave equation: 

2 

V 2 0 = h ^ f (2.02) 

0 TT* 
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The'scalar velocity potential is assumed to vanish as r-^oo and the 
Sommerfeld Radiation Condition is assumed to bo valid vrLth respect to the 
radial coordinate* 

The time dependence of the scalar velocity potential is taken to be 
= tf{r, e,z)eJ u,t as discussed in Section I* The scalar velocity 
potential $(r 9 e,z) must then satisfy the Helmholtz V/avo Equation 
(V^/k^) fS - 0* (2*02) An eigenfunction which is a solution to the 

Helmholtz Wave Equation, satisfies the vanishing boundary conditions and 
the Sommerfeld Radiation Condition, is developed in Appendix 1* 



This eigenfunction is stated: 

d ( \ - a sin(ne) =j2TT J j (2) 

K? ~ V? cos(ne) 



v Vd 



Where 



y« 2 a k 2 - (2 n f ) z 



(2.03) 

(2.04) 



/ = l|k 2 - (2 Trf ) 2 ’ 



(2 7rf) z < k 2 
(2 ff f) z > k 2 



^ = -j filr r ) 2 - k 2 " 

It is now proposed to determine a $ ( r,e,z ), constructed from a sum 
of the eigenfunctions (r,e,z), which will satisfy the boundary condi- 

tions at the surface of the cylindrical array* The technique to be 
employed in determining the complete scalar velocity function is a 
logical development from the technique of L* V* King used in the solution 
of the radiation impedance of a piston in an infinite rigid baffle* The 
constant A^^ transforms into *n( f ) as the eigenvalue f becomes a 
continuous variable* 



The sum of eigenfunctions is taken to be: 
^o© 
oo 

0(r,e,z) a 



n°0 



^ cos(ne) H n (2) (/r) e j Z7r T 2 df 



05) 



-oo 



°L. V. Kihg, On the Acoustic Radiation Field of the Pie j:6-Electric 
Oscillator and the Effect of Viscosity on Transmission, Can.Journ, Res., 11, 
pp 135-155, July-Dee .1934 
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The velocity in tho radial direction iss 
v r (r,ez) - - -gyr $ (r,e,z) 



oo 




n=0 




( 2 ) ( 2 ) 

f n (J>) Wl (A 7 ^rV r > cos(ne) .W* 

(2.06) 



The radial velocity distribution of the finite array will now be 
investigated# From the boundary conditions as defined, it is possible to 
write an expression for the radial velocity amplitude at the surface of the 
array: 



v r (a,e,z) 



oo 




n=0 



g(z) A n cos (n® ) 



(2.07) 



The funotion g(z) is equal to zero on the rigid baffle and on the i^* 1 

cylinder is equal to e^i. If the two expressions for the radial 

velocity, (2.06) and (2.07), are compared; the equality of the two 

quantities is possible if they are equal term by term, that is: 

( 2 ) ( 2 ) 

eUjAn. \ f n (J») * d/ 

( 2 . 00 ) 



-OO 



Define a new quantity. 



( 2 ) ( 2 ) 
^■/l(|Aa) ~ ^n-1 (,A a ) 
2 



A 



(2.09) 
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Then 



-\00 



g(z) 



g( Y ) o -/jZfry z d y 



( 2 . 10 ) 



-CD 



It can new be seen that G ( f ) and g(z) are a Fourier Transform Pair. 
Therefore : 






G( f ) « \ g(z) e"J 27r ^ z dz 



(2.11) 



-oo 



'b-yS 

(0) Q-i Zir $ z dz -f \ V 1 eJ^l ^ z dz 

-oo ^ "b-fZ 

(2i-l)b-/,S 

^ Vi e^i e-j ZTf t z dz / 

(2i-l)b-/Q^ 

p (2m-l )b-//3 






^ / 



N 00 



V m ej^ m e -J2rry 2 d2 ^ 



(0) e"j 25r ^ z dz 



(2m-l )b-yS 



( 2 m-l)b-/^ 



G(J) ; v x 1.-32"- ]T 1> -f 



i ....4 Vi e^i e-J^r < 2i -D b y 

V 

V V m e-J 2Tr ^ (2m-l)b 

vT 



m 



G(J) = 



i=l 



V. e“0 27r ^ (2i-l)b 



sin 2JT%£ > 

^jr 



6 



( 2 . 12 ) 



It is now established, after substitution: 

m 



?nW = 



2A n sin(Ziffp) V± e^ 1 e"j 2 ^ ? ( 2i “l) b 

i=l 



ru) — — no — ~ 

LVl^' a) " H n-l^ a )J 



(2.13) 



and 



/( r,e,z)s 



oo 



n=0 



oo 



X 



V* 

L-oo 



( 2 ) 

2A n cos(ne) siu(2^ )H n i^x) 

rt2i i2) v 
/“MT |Vi^ a> - Vi/ a >] 



m 

Vl e^i e "J2 (2i-l)b e j2/Tjf» z 

its 1 

(2.14) 

To obtain the average pressure, one uses the force equation 

vp /far : 0 

Sinoe V ■ - XJ ft and for the condition of time dependence e^ 1 ^, the 
force equation takes the form: 

p(r,e,z) « jw^(r,9,z) 

The average pressure amplitude on the active area of the h^ 1 cylinder is 
defined as p^. 

v/°c (2h-l)b+y3 

0(a,e,z) de dz 



Ph= J 



• Co 



4°tJ3 



-cx. (2h-l )b-yS 



(2.15) 
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The averaging process is carried out and the expression for the 



average pressure is, upon substitution for A n : 

oo 



oo 

Ph=j w ^° ^ / 

n»0 



n 

sin~^n«-) 

<n«5 



2 H 

n / 

^2 _ ( 2 ) 



(27rfp ) 



-OO 



oo 



V. e^e-J 4 ^ ^ ‘ (i “ h)b 



x y Vi O'" -*-6 



i-1 



d T 



(2.16) 

The average radiation impedance per unit area of the active area 



of tho h b ^ cylinder will 


be defined as 


< z > h* 


<Z> h = 


_ik_ 






V OYh 






h e 




Further define G~ « 


IT / 4 




<r z 


^1 - (*T ) 2 ’ 


(Vf ) 2 < 1 


<r z 


j- ^( Af ) 2 • 


• 1 (Vf ) 2 >1 



(2.17) 



(2.18) 

\ 

Then the radiation impedance per unit area of the active area of the 
th 



h cylinder may be written: 



oo 



< 2 >h=i 4 fi» 0 ^ 



n=0 



sin (n<x) 
'n p 

( nx ) 2 



P ( 2 ) 

sin (2 wj-^g ) (<r ka) 



-OO 



8 



X 



m 




V, 



e - 3 4>f (i-h)b d jo 



(2.19) 



The radiation impedance of an array of finite cylinders which pulsate 
uniformly with arbitrary amplitude and phase over a sector of each cylin- 
der has been determined* Evaluation of equation (2*19) involves numerical 
integration of an infinite integral* The integrand of (2*19) must be 

carefully examined in the region of -i-* A further handicap to 

A* 

evaluation of (2*19) is the fact that second order Hankel functions of 
complex arguments must be evaluated* These qualities of (2*19) are 
especially undesirable if tho radiation impedance is to be evaluated by 
digital computer techniques* Equation (2*19) will now be reduced to a 
form suitable for numerical integration* 

The following relation is stated and used in deriving a new form 
of equation (2.19)* 4 

K v ( jx) a - i-ir e "ji 7r ( v_1 ) h v ( 2 )(x) 

1 




oo 




n«0 






sm 






H- 



( 2 ) 



■n 



( CT^ka ) 



! , r (2) (2) 



m 




i = l 



v i „ j (V'i-lh) 



cos £ 4 It ^ b( i-h )J 







McLachlan, Bessel Functions for Engineers, Oxford at the Clarendon 
Press, 1955, p. 204. 
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oo 



oo 



■ti* ¥■ 2 Sn 

n=0 






i 

1_° 

-K 



sin^(2 Tk j* /3 )(-j) Kn(.i<Tka) x 

(2^(f/ 3 ) 2 [x^ (j<Tka)-/K n . 1 (jrka)] 



m 



X / -h- e J ) oos |- 47r ^ b(i-h)~| d f 



i=l 



( 2 , 20 ) 



The integrals will now be considered separately. Let 



1 



II = 



sin 2 ( 2 ^^) 4 2 ^(fka) 



0 



o r (2 ^ _ < 2 ) 

(2 irjp) 2 <r [Hn^i^ka) - H^^ka)] 



oos [4^^ (i-h)bj d 



(2.21) 



Let s sin £ 

then when J: 0 ; j = 0 

X=f; f-i w 



«" = - ( *■? ) 2 ' = oos f 

dj a COS ^ 



A. 



! 



Therefore ; 






sin 2 (k^3 sin | ) 

(k^sinj ) 2 Hn^(ka COS /) - Il<5}(ka COS /) 



^n^ 2 ^(ka cos J ) 



0 



cos 2kb ( i-h ) s in jj d ^ 

~ 



( 2 . 22 ) 
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Consider : 



oo 



. 2 



I 2 = 



1 

x 



Let 




in (Z'y . f/S) (-.1) K n (,i<rka) 

(Z 7 rf? )V [i^ (jr ka) T^n-lQ'*' ka)j 



cos [4irJ(i-h)bJ 



X'f B o SC ^ 



Then -when - 1_ ; ^ a -g-H' 

X 

T - OO ; / 5 0 



<r Z -3 |(XJf ) 2 - 7 1 s -j ootf 



dj = -esc/ cot f 



•< <LL 

X 



sin^(k ^ esc / ) 



K n (ka c ot / ) 



(kyO ) 2 cso/ (ka cot f ) K n-1 (ka cot/ ) 



(2.23) 



cos ^ 2kb (i-h)c sc/J d^ 



(2.24) 

The expressions (2.22) and (2.24) are substituted in (2.20) and the 
resulting equation for radiation impedance per unit area z ^ takes the 
more tractable form: 



-!r 



oo 



<z> h 3 J d (°o G J? 




n=0 



sin^(noc ) \ sin^(k/sin/) n 



H ^)(ka cos / ) 
n 0 



(n«,)‘ 



(k/,SinS>) ^-o S p-.l<!/ k a oosj) 
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m 



cos [2kb(i-h)sinj]dj 



i=l 



oo 



/ j Q (c 



n=0 



h 

V 

2 






c K ,i sin 2 (not.) 



sin^ (k/3 e sc ^ ) Kn(ka cot ^ ) 

(k^6 ) esc f K^-^ka oot£) /K n- ^(ka cotj) 



m 




i=l 



(2*25) 



The evaluation of <z> h ^ as b 0en reducod to that of evaluating finite 
integrals and functions of positive real variables for which mathematical 
tables are readily available# A further manipulation of (2*25) will be 
made in order to separate the real and imaginary components of < z> ^# 

Let 



B n (ka cos ^ ) e“j^n(ka cos £ ) 



II n ( 2 )(ka cos / ) 



II 



( 2 ) , ( 2 ) 

y^(ka cos ? ) - H n _^(ka cos^ ) 

(2.26) 



Then the expression for < z> ^ transforms into the following expression: 



12 




^n-/l (ka cot > )-/K n «.i(ka cot^ 



The radiation impedance per unit area of an array of m finite 
cylinders which pulsate uniformly over a given sector of each cylinder 
with phase and amplitude arbitrary but constant for a given cylinder is 
expressed by equation (2.27). There are no restrictions stated on the 
wavelength of the radiated sound, on the cylinder length or radius, or 
on the cylinder spacing# Although equation (2.27) is for a system with 
uniform velocity and phase across each radiating element, the method 
of solution is not limited to such a velocity distribution. Any 
velocity distribution in azimuth which can be expressed in a Fourier 
Series is possible of solution using such technique. Further any 
velocity distribution in the z coordinate for which a Fourier Transform 
exists is also capable of yielding a solution. The distribution in the 
z coordinate need not be tho same on each cylinder for a solution could 
be constructed piece by piece for such a system. 



14 



3. Spooial Cases 



The validit/ of equation (2, 27) will be considered by investigating 
two special cases of the finite array for which solutions already exist. 
The first array to be considered is the single uniformly pulsating 
cylinder in an infinitely long, rigid baffle* Let the array be con- 
structed as follows: 

1* There is one active cylinder of length - 2 f> and radius = a. 

2. The cylinder is mounted coaxially in an infinitoly long, rigid, 
cylindrical baffle of radius = a. 

3. The cylinder pulsates uniformly with radial velocity = 



The Fourier Series equivalent of the third condition as stated above 
is expressable in one term. That is, on the active cylinder tho radial 
velocity may be expressed as one term of the Fourier Series (2.0l). 



With reference to equation (2.27), tho above statement or condition 
reduces the summation with respect to n to one term, namely n = 0. The 
condition of one cylinder reduces the summation with respect to i to one 
term, namely i = 0© Let the radiation impedance per unit area of the 
uniformly pulsating cylinder be < z> • Then from equation (2.27), 



V 0 








0 
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JT 



\ 



2 



■/ ^ 



0 




(3.01) 



Equation (3*01) is written with substitutions according to (2*26): 



radiation impedance of a single uniformly pulsating cylinder in an 
infinitely long, rigid cylindrical baffle if the variables of integration 
are changed and it is allowed that Robey’s equation contains an error in 
sign# (This is demonstrated in Appendix IX)# 

It will now be shown by uso of the Dirac Delta Function that 
equation (3*02) will transform into the well known expression for the 
radiation impedance of an infinitely long, uniformly pulsating cylinder# 
Equation (3*02) is written in the form of equation (2#19) by trans- 
formation of variables: 

H# Robey, On the Radiation Impedance of Finite Cylinders, J. Ac oust. 
Soc. Am#, 27, pp. 706*711, 1955 





T 

2 



0 



(3#02) 



Equation (3.02) compares with Robey's equation (18) for the 
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^oo 



< 2 > = 3f 0 ° 



sin 2 (2 7T/3 f ) Ho^V^ka) 



-oo 



^{Zrr/3 )yV„ (2) 



d r 



H, (<Tka) 



(3.03) 



Consider the expression 

a sin 2 ( 27 7 d> 

2»y»*oo 2T//3 ) 2 



(3.04) 



Expression (3.04) has the following qualities: 

1# It is zero for all j° except J 3 - 0. 

2. At if - 0, (3o04) tends to infinity. 

3. The integral of (3.04) from f equals minus infinity to 
^ equals plus infinity is equal to unity. 

The above is tho definition of a Dirac Delta Function. Therefore: 

■ 2 , 



lim 

2lTp->00 



-L. ? in . ( 2 r A 1 r } = 8 (f) 

TT '1 , „ 4 o <> 



2 */3 ( f Y 



(3.05) 



and 



oo 



lim <z> 

2irj3-» , c>o 



= ^o c 



^ cr ) 



(erka) 



<r h- 



( 2 ) 



(<T ka) 



d r 

J 



-oo 



H 



( 2 ) 



(ka) 



lim <z> = j£ Q 0 -r-r 

■A— k no .. 



2lT|'i~*- OO 



H, 



(ka) 



(3.07) 
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The equation (3.07) which is the radiation impedance per unit area 
of an infinitely long, uniformly pulsating cylinder compares with the 
well known expression for such radiation impedance.® This concludes the 
discussion of the validity of equation (2.27). 

£ 

L# Beranek, Acoustic Measurements, John Wiley and Sons, Inc#, 1949, 

p 0 60 
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4. Beam Steering 



The application of equation (2# 27) to 5 specific array will new be 
considered in order to study the effect of electrical bean steering on the 
radiation impedance of the array. The array is to be constructed as follows 



kb = 0 

2. The cylinders are mounted coaxially in an infinitely long, 
rigid cylindrical baffle with lea = Y 

3. The cylinders pulsato uniformly over their surface (symmetry 
in azimuth). 

It is proposed to investigate the effect of electrical beam steering on 
the radiation impedance of the individual cylinders. Electrical beam 
steering will be accomplished by phase leading each cylinder a uniform 
number of degrees ( i ) with respect to its position in the array. 

Equation (2.27) is written to conform to the above boundary 
conditionso 



1. There aro nine active cylinders with kp - -g-TT', ka and 



r 

02 




9 



< 2 >h = j 



0 



-/ j 




0 
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(4.01) 



The result of the evaluation of (4*01 ) by numerical integration 
using the trapezoidal rule and A f = — - is plotted in figures 2 and 3 
for cylinders number one, three, and five# 

For such an array, the electrical beam steering has but a moderate 
effect on the radiation resistance of any element. Electrical beam 
steering affects the radiation resistance of tho end elements of tho array 
to a greater degree than it affects that of the center elements. 

Electrical beam steering has a more drastic effect on the radiation 
reactance of such an array# Radiation reactance in this array varies as 
much as ten to one for the center element of the array# 
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APPENDIX 1 



The Helmholtz Wave Equation is restated in cylindrical coordinates: 

/ -i 4 h + ¥4 v V : ° < ai - oi > 

d r 3 r d o 



The solution to the above equation is assumed to exist in a variables- 
separable form: 

0(r,e,z) = R(r)G(e)Z(z) (A1.02) 

(Al#02) is substituted into (A1*01) and the resulting equation is: 

1 d^R 
R dr‘ 



/ 1_ _dR_ / 1_ d 2 fr , 1 d 2 Z ,2 

*2 Rr dr ftr2 d« 2 Z dz^ 



(A1.03) 



1 d^Z 

It can now be seen that — — * is equal to a quantity which is independent 

Z dz^ 

of z; therefore the following is true: 



-l — sL«£ - - (2 ■fr $ ) 2 



(A1.04) 



Z dz2 

Where Jp is an eigenvalue to be determined by the boundary conditions 
assigned to the array* The solution for Z(z) may now be written# 



Z(z) = 2 /by Q z 

Continuing with the solution, it can be seen that 

to a quantity that is independent of e; therefore: 

-L = . n 2 

9 d9 



d 2 9 



9 de 



(A1.05) 
is equal 



(A1.06) 

Where n is an eigenvalue to be determined by the boundary conditions 
assigned to the array* It will be found more convenient to write the 
solution for 9(e) as follows: 

0(e) = c n cos(ne) / sin(ne) (A1.07) 
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The assignment of arbitrary eigenvalues to the solutions for Z(z) and 9(e) 
determine the form of solution which ( r ) must assume# 

d \fy 1 <jjk 2 _ (21 r'S ) 2 ] n 2 j, _0 (A1.08) 

dr r dr r 

In the following discussion, the quantity k 2 - (2ft^ ) 2 will be defined as 
equal tO/^~^« The quantity /^ must reduce to k as (2^^ ) 2 approaches zero 
therefore : 

= ^k 2 - (2 ~7r'f ) 2 ' (2 TTf) z < k 2 

-j J(27r^ ) 2 -k 2 (27?J) 2 > k 2 (A1.09) 

To satisfy the Sommerfeld Condition of Radiation, the solution for R n ^ ( r ) 
must be : 

R n> (r) = H n ( 2 ) <yur) (A1.10) 

The expression for R n ^ (r) as r approaches infinity is: 

R ( r ) r 1 2 e _j( / r (Al.ll ) 

This expression coupled with the time dependence eS 14 ^ represents an 



out g o ing wave # 

The wave equation has been solved consistent with the Sommerfeld 
Condition of Radiation# It is now possible to state a eigenfunction which 
is a solution to the wave equation; and if the boundary conditions imposed 
on the array are consistent, a sum of the eigenfunctions should satisfy 
the additional boundary conditions# 

The eigenfunction is stated: 



a* <*•••*>■*»* 



(A1.12) 
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APPENDIX II 



To substantiate the statement that equation (18) of Robey's paper is in- 
correct, the development of equation (18) from equation (16) will be inves- 
tigated as yrell as the validity of equation (16). Further it will be shown 
that equations (19) and (20) of Robey's paper are correct. 



z qq = 



,oo 



H 0 (2) (Ka) sin 2 (jL) 
K H-^^Ka) 



dj 



(16) 



Where K = (k 2 -f 2)« 
Rewrite equation (16) as follows 

7^p° 

4 a feu) 



4^L 



sin 2 (/°L) 



H, 



( 2 ) 



(Ka) 



(ff K Hp( 2 ) (Ka) 



df (A2.01) 



-oo 



Where < z^q> = ‘the average radiation impedance per unit area 
Consider the following expression: 



lira BAn - 2 IL L i 
L-^oo ff L( j° ) 2 



(A2.02) 



The following facts are true: 

1. (A2.02) approaches infinity as approaches 0. 

2. (A2.02) is zero for all other f . 

3. The integral of (A2.02) from J 5 = -oo to 6P a -/oo is equal to unity. 
Therefore (A2.02) is a Dirac Delta Function, and (A2.01) may be written 

as follows: 



lim < z qq > * 



L~>oo 




i ( f) H " (2)<Ka) 



K 1^(2) (Ka) 






(A2.03) 



-oo 
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Thus the radiation impedance per unit area of tho uniformly 
pulsating cylinder as derived by Robey reduces to the following as the 
length of the cylinder approaches infinity: 



iim Cz qq > = -jfo° 
L-* oo 



Hq ( 2) (ka) 
Hj < 2) (ka) 



(A2.04) 



This does not agree with the well known result for an infinite 

II Q ( 2 )(ka) . 



cylinder* Further the angle associated with " Q "/ ' is always in the 

III (^) 

fourth quadrant; thus the radiation resistance given by (A2#04) would 
always be negative# 

It oan be assumed that Robey’s equation (16) is correot except for 
sign, and it will be shown that Robey’s equation (18) has an error in 
sign but that equations (19) and (20) are oorrect# Equation (16) is 
restated with the correct sign: 



z qq 




sin^ (J°L) H c ( ^ (Ka) 
</> 2 ™i < 2 > ( Ka ) 



(A2.05) 



z qq = d 8a ^ 



J 



sin 2 (f L) Hq ( 2) ( f k 2 -y 2 ’ a) 

(1 ) 2 |k 2 -f 2 ' Ilj^H |k 2 - f S a) 



4 




k 



sin 2 ( J° L) 
</) 2 



(-j) K 0 ( \ 


J 2 -k 2 ' a) 


(-j)\|/ 2 -l 


^ K l( ^ / 2 -k 2 a) 



dy 



(A2.07) 
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To reduce (A2.07) to tho form of Robey's equation (18 ) consider the 



following integrals: 
k 

*1 = 



S in 2 (y L) H o (2) ( I**- X 



(f ) 



Hl {2, ( [7- J> 2 'a) 



dj> 



Let 'f = k sin T?' 
then when = 0, "tf" s 0 

Jn k, -XT - nf 
2 



and |k 2 ^ C os^ 

dy a k cos - ^ dz*' 

Therefore 



'2 






sin^(kL sinz?) H 0 ^^(ka cost?") 



(siniy)' 



(ka cos t?) 



d & 



Consider 



*2= 



oo 



sin 2 (J > L) 



U 

k 



(jr 



K 0 ((y 2 -k 2 a) 



2 -k 2 a) 



Let y : k cosh/ 7 

then when If : k, « 0 
y = oo, oo 

and U z -k^ = k sinh^ 

df = k sinh^ d<y 



df 



(A2.07) 



(A2.08) 



(A2.09) 
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Therefore 



o© 



X 2 = ~*2 



sin 2 (kL cosh ) K p (ka sinh ) 

(cosh-/') 2 K x (ka sinh ) 



d f 



(A2.10) 



Equation (A2.05) may now be written as 

df 
'2 



z qq = J 



8 a (j-) 



r 

j 



o 

/- o 






8a 'P n cu 



) 



2 (2) 

sin (kL sini^) 1I Q y (ka cos ) 

(sinx *) 2 oos^) 



sin 2 (kL coshy/) K p (ka sinh^' ) 
(cosh/' ) 2 K-^ka sinh/-) 



■/ 



d'f (A2.ll ) 



Thus from (A2.ll), it can be seen that the terms in Robey's equation 
(18) should have identical signs; and from the initial argument about 
equation (16),* the terms should be positive. 

Robey's expressions for $ and & as they appear in equations (19) and 
(20) will now be derived from (A2.ll). The C/ torn is self explanatory and 



correct. Examine the expression 

H n ( 2 )(ka oos "O' ) a H q ( 2 ) (x) 



II- 



•j_ ^ 2 ^ (ka cos ~iy ) 



V*' (x) 



( 2 ) 

J 0 (x) - jYo(x) 
J].( x ) - JY 1 (x) 



where x = ka oos 



Jl(x) 4 JYlU) 



JqU) / Y]_(x) 

J Q (x)J 1 (x) 7 ^Y 0 (x)Y 1 (x) ^ J 0 (x)Yi(x) - J 1 (x)Y 0 (x) 

J 2 (x) /Y 2 (x) J 1 (x)/Y 1 (x) 
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J 0 (ka Gos^)Ji(ka cos /Y 0 (ka cos^) Y\ (ka cos <^) 




J 



. 2 



ka cos 




Thus it can be seen that Robey 1 s coefficients Jc , $ , and £ in 
equations (19) and (20) are correct as stated but are not consistent with 
equation (18)o 

The above arguments are believed to be sufficient to show Robey 1 s 
equation (18) is not correct as stated, but would be correct in the form 
of equation (A2.ll)* Equation (A2*ll), if stated in terms of average 
radiation impedance and with a change of variable, is equivalent to 
equation (3#02)# 

7 

Go N. Watson, A Treatise on the Theory of Bessel Functions, The MacLiil 
Company, 1945, p 77 
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(Figure J. ) 
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